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Abstract 

Wireless Body Sensor Nodes (WBSN) are frequently used for real time IoT-based health monitoring of 
patients outside the hospital environment. These WBSNs involve bio-sensors to capture signals from a 
patient’s body and wireless transmitters to transmit the collected signals to a server located in private/public 
cloud in real time. These WBSNs include hardware for processing of signals before being transmitted to the 
cloud. Simultaneous occurrence of all these processes inside energy constrained WBSNs results in considerable 
amount of power consumption, thus limiting their operational lifetime. Due to the inherent error-resilience 
in signal processing algorithms, most of these data reaching the servers are redundant in nature and hence 
of not much clinical importance. Transmission and storage of these excess data result in inefficient usages of 
transmission bandwidth and storage capabilities. In this paper, we develop a real time encoding scheme that 
performs iterative thresholding and approximation of wavelet coefficients for sparse encoding of bio-signals 
(ECG signals ), thereby reducing the energy and bandwidth consumption of the WBSN. The encoding scheme 
compresses bio-signals (ECG signals), while still maintaining the clinically important features. We optimize 
various process parameters to model a low power hardware prototype for the implementation of our algorithm 
on a real time microcontroller based IoT platform that operates as an end-to-end WBSN system in real time. 
Experimental results show a system-level energy improvement of 96% with a negligible impact on signal 
quality (2%). 

Keywords: Wireless Body Sensor Nodes; IoT-based Health Monitoring; Low Power Hardware Prototype; 
ECG Signal; Discrete Wavelet Transform; Sparse Encoding; Approximate Computing 


1. Introduction 

In the present world, the growth of cardiovascular diseases and interests in personal healthcare fuel up 
the requirement for a system that provides accurate and non-stop medical supervision mm- For chronic 
cardiovascular patients, the traditional healthcare systems entail two categories of monitoring, which are i) 
holter monitoring that requires a smaller time range of about two days, and ii) event monitoring that spans 
across a week. These monitoring systems involve bulky and obtrusive sensor devices like Electrocardiogram 
(ECG) sensors that limit patient autonomy [3]. Moreover, recording over multiple channels and in multiple 
sessions with sampling rate varying from 125-500 Hz, leads to a huge amount of data generation [4j. Thus 
continuous medical supervision is unsustainable for traditional healthcare system due to increasing costs 
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resulting from energy consumption by the sensor nodes (devices), data transmission, and larger storage 
requirement [5]. 

Wireless Body Sensor Nodes (WBSNs) have emerged as a feasible/portable solution to enable the long 
term acquisition of bio-signals, outside of a hospital environment mm- WBSNs are embedded sensor appli¬ 
ances consisting of low power nodes aiding in the acquisition of ECG signals Bang. Smart WBSNs incorporate 
various innovative signal processing applications to reduce the ECG data to be transmitted HD). The am¬ 
bulatory ECG data collected is affected by noise like powerline interference, electrode motion, muscle noise, 
baseline drift, etc HD. WBSNs, operating on battery, are usually heavily energy-constrained that requires 
the development of energy and storage-efficient embedded signal processing and data compression algorithms. 
These processing algorithms should therefore be simple and of low complexity [12]. Also, the overall system 
latency required to process and send data should be minimal as delay cannot be tolerated at the time of 
crisis mm- 
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Figure 1: Wearable body sensor node architecture 

A typical WBSN comprises the following three subsystems depicted in Fig jT| i pinp ): 

• Sensing Subsystem: This comprises biosensors and analog to digital (ADC) converters. Analog biosig¬ 
nals acquired through the biosensors are converted to digital form for processing in the next subsystem. 

• Processing Subsystem: Noise filters and signal processing units with embedded algorithms for signal 
compression and energy optimization, form this block. Digital biosignals from the sensing subsystem 
are processed and prepared for transmission by the next subsystem. 

• Communication Subsystem: This comprises the Ultra Low-power WiFi module which is used to transmit 
the processed signal data to a medical server. 

As demonstrated in Fig. [T] the processed digital signal, along with the other relevant details of a patient 
are transmitted to a medical server. Finally, this server sends the patient’s data to the doctor or healthcare 
provider or any other authority under concern ma¬ 
in order to perform the desired processing, the ensuing smart WBSNs embed low complexity computing 
capabilities that have non-negligible energy footprints HU. Note that the design considerations for WBSNs 
considered in this paper are a) power consumption, which includes both transmission power and execution 
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power, b) portability, c) real-time availability, d) energy efficiency and e) scalability. Energy consumption of 
the system depends upon the sum of the energy of the individual components in the hardware platform such 
as the sensing, processing, and communication subsystems as shown in Fig. |T] However, the energy efficiency 
can be examined as the energy bottleneck through embedded digital signal processing (DSP) algorithms [15], 
Real-time availability depends upon the time complexity of the DSP algorithms and the hardware latencies 
jBj. In the case of healthcare application, a real-time system is actually a soft real-time system in which 
some latency is permitted [19 . Increasing the energy efficiency and performance are the two major driving 
factors considered while designing a real-time IoT based health monitoring system. For energy-tolerant 
application such as ECG signal analysis, guaranteeing 100% correctness is not always mandatory and we 
leverage this error-resilience to gain system-level energy significantly in this work. In order to achieve this 
gain in the system level energy, this work first investigates the error tolerances of our proposed algorithms 1 
and 2 (refer to Section [4]) through regulation/optimization of certain design parameters that are subsequently 
used to formulate an energy saving methodology using approximate computing. The experimental results 
demonstrate the significant system-level energy gains that can be achieved with negligible impact on the 
clinical value of the transmitted signal. 

A WBSN based device continuously monitors ECG signals from a patient’s body. This continuous moni¬ 
toring creates excessive volumes of data that, in turn, requires larger storage space and incurs higher trans¬ 
mission cost and power consumption. Also, this large amount of data requires extra processing time to be 
executed by the Processing Subsystem of the WBSN device, resulting into massive power consumption. Thus, 
frequently data compression algorithms are utilized for sparse encoding of this data. As a result of sparse 
encoding, the processing requirements decrease, which result into less power dissipation in the Processing 
subsystem. As the number of transmitted bits (that forms the packets) is decreased, the transmission is 
decreased too. At the same time, since sparse encoding algorithms are lossy in nature, the signal recon¬ 
struction at the receiving end will not be exact, and will result into some deviation between the original 
and reconstructed signals. However, in clinical domain, maintaining a very low percentage of root mean 
square deviation (PRD) among the original and reconstructed signals, is a matter of utmost importance. For 
example in this case, let x be the original ECG signal and x re con be the compressed version of x, then PRD 
is a measure of deviation of x recon from x and is given by Equation 1. If PRD will be high, certain subtle 
clinically important features in the signal will be lost; this will be inappropriate in this case. Increasing the 
sparsity of the encoded signal would increase the PRD and in in turn would lessen the power consumption, 
storage requirement, transmission time and system’s computing overhead. The sparsity achieved by the re¬ 
constructed signal is measured by the Compression Ratio (CR) which defined as the ratio of number of bits 
required to encode the original signal ( B a ) to the number of bits required to encode the compressed signal 
(. B c ) as described in Equation-2. If the reconstructed signal becomes too erroneous, it will have a high PRD 
and hence it will lose its clinical importance to some extent. Although, it is important to reduce the PRD, 
occasionally with a very small increase in it, we can get significant levels of energy saving. Sometimes a 


3 


very small PRD may not be even required to obtain correct ECG analysis, and in those cases the PRD can 
be traded-off for much higher amounts of system-level energy savings because of improved CR. By solving 
an optimization problem, this paper aims at establishing a trade-off between the PRD and the compression 
ratio/energy consumed. 


PRD{%) = l|X x 10Q 

Imh 

(i) 

Cfl =t 

(2) 


To explain our idea, we use illustrations using real experimental data as shown in Fig. [2] that shows both 
microcontroller execution time ( mctime ) (energy consumption of the processing subsystem) and transmission 
time ( txume ) (that impacts the communication subsystem energy energy) decrease together with an increase 
in error bound (e r ) ( PRD is obtained by representing (e r ) in %) between the original and reconstructed 
signals. The value of e r is increased from 0.0 to 0.6 that is depicted in Fig. §> to Fig. W- With an 
increase in e r , the sparsity of the encoded signal increases, as a result of which mc t i me decreases. As e r 
increases, few numbers of bits (or packets) are transmitted and hence, txtime decreases. Note that mctime 
and txtime are directly related to the microcontroller execution energy ( mc energy ) and transmission energy 
( tx ener gy ), respectively. However, as we see, both mc energy and tx energy are indirectly proportional to e r . 
In Fig. i- and Fig. ih we see that as e r increases by 0.2, both mc t ime and tx t i me decrease by 41 ms and 
1.1 ms respectively. Decrease in both mctime and txtime lead to an overall decrease in energy consumed in 
both processing and communication subsystems. This in turn results in a reduction in total system energy 
given by tot energy =mc energy +tx energy . Note that the sensor energy ( sen energy ) is assumed constant and also 
turned out to be very small as seen from our energy consumption data and hence, has not been mentioned 
explicitly in the total energy calculation for this paper. Thus, from Fig. - Fig. [2]T, a tradeoff between the 
amount of tolerable error and the energy consumption can be observed. To achieve this system level energy 
improvement with negligible impact of signal quality, the following algorithms have been proposed in this 
paper: 

• Encoding Algorithm: that can directly control the sparsity level and hence, the amount of error intro¬ 
duced in the reconstructed signal with respect to the original signal. 

• Optimization Heuristic: that aims to minimize the the total amount of energy consumption of the 
WBSN while incurring negligble impact on application-level quality. The tolerable small quality/PRD 
degradation is stated as an input design parameter. 

This paper presents one example of IoT based ECG system that we use to demonstrate how approximate 
computing can be used to reduce the energy consumption drastically of these types of IoT-based medical 
systems. There will be available other similar systems too, and our generic approach is to be applied to those 
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Figure 2: Illustrative Figures (a)-(d): Variation of total time and power consumption for various values of introduced errors. 
Error increases with an increase of sparsity from (a) to (d). 
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systems as the approximations are performed both at the system (hardware) and the application (software) 
levels. As stated earlier, although it is mandatory to reduce PRD, occasionally with a very small increase in 
it, we can get significant levels of energy savings. Sometimes, very small PRD may not be even required to get 
correct ECG analysis, and in those cases the PRD can be traded-off for much higher amounts of system-level 
energy savings as performed in a systematic way in this paper. We introduce our proposed error- resilient 
algorithms for demonstrating the usefulness of the approximation which would lead to energy saving in the 
designed system. These values can be set initially based on inputs provided by medical and signal processing 
experts and then these values can be adapted accordingly based on the ECG signal characteristics and the 
specified error bounds given by the clinical experts. 

2. Related Work and Our Contributions 

With the increase in health awareness among people, new demand for real time remote health monitoring 
had been observed recently mm- Information and Communications Technology (ICT) had been able to 
cater for this demand significantly [22]. Over the years, ICT improved the efficiency of healthcare sector 
at a great extent [23, [24]. This led to the growth of e-health devices: devices which were used in ICT for 
autonomous monitoring and controlling of healthcare related issues [211 OH EH EH] .However, there are certain 
concerns related to the application of ICT in healthcare service Gum- The major issue is that e-health 
devices consumed enormous energy, much of which was wasted El HUE]. These devices were connected 
with a large number of sensors and hardware components for information processing and signal transmission 
[521 [551 . Real time health monitoring demanded these e-health systems remained to be powered on, to collect 
data in a continuous manner, for longer periods of time. However, not all collected data are clinically useful, 
and only a fraction of the information collected was utilized by the concerned authorities for making decisions 
[54] . Thus, the wastage of energy was enormous in these devices [35) 36]. Another issue with these devices 
was to mention that they generated a high volume of data which often required larger storage [53 [3S]. This 
demanded for efficient real time compression scheme mm- While most of the existing worked on signal 
compression either achieved low PRD between the original and reconstructed signals or high compression 
ratio (CR) [ill 1521 H51 13H 05] .Compression Ratio is the ratio of number of bits required to encode the original 
signal ( B a ) to the number of bits required to encode the compressed signal ( B c ). But, the main problem lay 
in achieving both low PRD and high CR simultaneously. Some previous works on ECG signal compression 
assumed arbitrary threshold values for sparse encoding of signals [461 1171 08] 09]. The real life situations 
showed that the changes of threshold values of ECG were depending upon heart conditions of persons under 
clinical observations. Thus, hard thresholding might result into loss of clinically important information [50] , 

In this work, we design a real time IoT based health monitoring system with optimized process parameters. 
The optimization heuristic provides solution scheme which allows to reduce power consumption, and improve 
data compression while maintaining significantly low PRD. In this paper, we propose a sparse encoding 
algorithm consisting of two schemes namely geometry based method (GBM) and the wavelet transform 
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based iterative thresholding (WTIT). The sub algorithm GBM reduces the minimal ECG voltage values to 
zero level. Subsequently, WTIT encodes the ECG signal in time-frequency domain, obtaining a high sparsity 
level. Our proposed encoding algorithm has the lowest PRD and achieves the highest CR, when compared to 
the previous ECG compression techniques that utilised compressed sensing .This comparison is demonstrated 
in Table |Tj The absolute time complexity of the whole procedure exclusive of the transmission media time is 
O(n), while most the previous related works had time complexities in the order of 0(n 2 ). 


Table 1: Comparison of our proposed algorithm with previous algorithms 


Algorithm 

PRD (%) 

CR 

Simultaneous OMP i5T] 

2.57 

7.23 

Basis Pursuit (1) [52j 

1.66 

4 

Basis Pursuit(2) [55] 

9 

3.44 

Basis Pursuit (3) (54j 

9 

2.5 

Our Encoding Algorithm 

2.27 

9.4 


3. Theory 




WBAN 

TIER 1 


Gateway Devices 

TIER 2 


End-user Healthcare Monitoring Application 

TIER 3 


Figure 3: IoT-based health monitoring system 


In the three tiered architecture of an IoT-based health monitoring system illustrated in Fig. [3j Wireless 
Body Area Network (WBAN) forms the first tier IoT platform of WBSNs. WBAN comprises miniature 
biosensors called WBSNs that can be both wearable and implantable. These biosensors collect physiological 
signal, in this context of ECG signal, and perform sampling, filtering and processing operations on the 
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Figure 4: WBSN circuit diagram 


data. WBSNs used in this work are smart sensors containing user modifiable section of signal processing. 
All the biosensors connected through short range communication form an ad-hoc sensor network. The 
second tier is constituted of the gateway devices like smartphones and laptops. The function of the gateway 
system is to provide the connection between ad-hoc sensor networks and infrastructure based Wide Area 
Networks (WANs). Ultra Low-power WiFi is used in this work to launch the transmission of the processed 
and compressed data from the WBAN to the Gateway Devices. It transmits information to one or more 
heterogeneous networks depending upon the application. The third tier comprises WAN and the number 
of end-user healthcare application providers as illustrated in Fig. [3j The end-user healthcare application 
provider has a processing unit and an user interface namely graphical user interface (GUI). Vital signs 
appearing in GUI are monitored real-time by healthcare professionals to keep track the health status of a 
patient under observation. Note that these WBSNs are highly energy constrained as they are embedded or 
implanted on or within the body of the human being. 

The IoT-based health monitoring system shown in Fig. [3| constructs Tier 1 nodes in WBAN that 
is segregated into three subsystems: i) Sensing, ii) Processing and iii) Communication as shown already in 
Fig- [TJ The power subsystem is abstracted in this system. The processing and the communication subsystems 
are the areas in which we develop energy efficient techniques in this work. The digital signal processing stage 
is user modifiable, and subsequently an inexact signal processing is implemented in it for achieving energy 
efficiency. In this work, we use a commercial-off-the-shelf (COTS) microcontroller for designing the sensing, 
processing and communication subsystems. The circuit design of a typical COTS microcontroller, namely 
WBSN in this paper, is shown in Fig. [4j 

• The Sensing subsystem comprises AD8232 ECG monitor with a lead setup of three electrodes. It is 
a dedicated front end integrated circuit heart rate monitor. It is designed to extract, amplify, filter 
small biosignals in the presence of noisy conditions, such as those created by motion artifacts or remote 
electrode placement. The monitor contains an ultra low power Analog to Digital Converter (ADC). 
The output signal is acquired on an embedded microcontroller. Input signals are taken from Right 
Arm, Left Arm and Right Leg forming an Einthoven triangle. 






































• To develop the Processing subsystem , we use Arduino UNO R3, which removes noise from acquired 
biosignals, adjusts baselines, encodes signals using sparse coefficients, and optimizes system parameters 
to consume minimum energy. Arduino UNO comprises a low power unit, an efficient microcontroller 
with built-in ADC and an USART communication feature. 

• The WiFi module in Arduino transmits real time data to a server located typically in public/private 
cloud. Those data will be retrieved by physicians or hospitals’ personnel. This Arduino based WiFi 
shield, ESP8266, is used in this system to design the Communication subsystem. 

While considering the DSP stage of WBSN, a practical scenario of embedded bio-signal analysis is assumed 
here for our case study. Embedded bio-signal analysis is a domain in which energy efficiency is the key, while 
the applications are inherently error-tolerant to a certain degree. In this work, we utilize this intrinsic error 
resilience of algorithm to achieve unprecedented level of saving of the energy of the system. We develop 
the algorithm in a way that it requires minimal hardware and a negligible energy overhead with real-time 
transmission. The bio-signal, considered in this work, is an electrocardiogram (or ECG signal). In ECG 
signal, the deterministic and stochastic feature are exhibited, and the repertoire of the methods applied 
for its analysis include (i) the extractions of morphological features, (ii) the analysis in time and frequency 
domain, (iii) the statistical methods and classification procedures. In order to study them, we develop, in this 
work, the parameter extraction method namely geometry based method (GBM) [55j , and wavelet transform 
based iterative thresholding (WTIT) (SB,. Utilizing the quasi-periodic nature of ECG signal, GBM increases 
the sparsity of the signal in the time domain. Due to localization properties in time and frequency domains, 
WTIT further approximates the signal in wavelet domain with a few coefficients creating a sparser signal, 
which is finally transmitted after having undergone Huffman encoding. Therefore, the proposed framework 
is divided into three sub-algorithms: 1) Noise filtering, 2) GBM, and 3) WTIT. 

3.0.1. Data Acquisition Through Noise Filters 

The processing section of WBSN begins with the wavelet filtering of ambulatory sampled ECG signals. 
For proper ECG signal acquisition, the first step is noise suppression performed by noise filters, and then 
the enhancement of ECG signal amplitudes through an amplifier circuit built of operational amplifiers. The 
filtering techniques used for denoising ECG signal involve linear and non-linear methods; either model based 
or model free approach. In this paper, noise filters employ wavelet based denoising (model-free thresholding 
method), brought it to completion by decomposing the signal using band pass filters corresponding to a 
mother wavelet basis function. The mother wavelet chosen should be maximally compatible with the ECG 
signal structure. 

Wavelet based denoising involves the use of Donoho’s universal thresholding. We set decomposition level 
dependent threshold value X) per Donoho’s universal thresholding m-- 

7) = log e l (3) 
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where, a is noise variance of the signal at each level of wavelet decomposing and l is the length of the 
analyzing signal. For a given signal u[n], let there exist an orthogonal matrix W such that the discrete 
wavelet transformed matrix U[n] is given by: 

U = Wu (4) 

Let, the additive noise to signal x[n ] be y[n], which has noise variance a 2 . Hence, the resultant signal z[n] 
is given by: 

z[n] = x[n] + y[n] (5) 

The inverse wavelet transform is computed as: 

u = W T U (6) 

The wavelet transformed matrix U[n] is modified to U[n\ by applying thresholding and as a result u[n] 
becomes u[n] 

u = W T U (7) 

Reduction of wavelet transformed coefficients produces significantly noise free estimate, which is done either 
by hard or soft thresholding. 

T hard {U-T l ) = UO{\U\>T l ) (8) 

T soft {U ;T,) = sgn(U)(U - T,)0( \U\ > T t ) (9) 

O is the usual indicator function and sgn represents signum function. In order to eliminate noise, we 
consider soft thresholding of wavelet coefficients because hard thresholding results loss of information from 
the original signal. The noise free ECG signal is free from zero mean band limited Gaussian noise i.e., the 
minute interference presents in the signal but it is not entirely free from the periodic baseline drifts. The 
noiseless ECG signal is processed using our proposed encoding algorithm. 


R R 



3.0.2. Geometry Based Method (GBM) 

Clinical assessment of ECG (Fig. [5]) relies on the evaluation of its morphological features in time domain, 
those features are position, duration, amplitude and slope of its QRS complex and segment. The first step 
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is the detection of the QRS complex which serves as a marker for averaging of heart cycles. To determine 
feature extraction from ECG beat, different approaches using Hermite Functions, wavelet transforms etc., 
were proposed. In this paper, we develop a time-domain based method GBM for the same purpose. GBM 
belongs to the preprocessing stage of Sensing subsystem that is shown in Fig. 1. This pre-processing technique 
is only applicable when ECG signal shows likely a periodic variation. In an ECG signal shown in Fig. s only 
P, Q, R, S and T waveforms and their relative separations are clinically important features. From preset 
threshold values and the rate of change of the amplitudes of these waveforms, GBM automatically identifies 
them. It then diminishes all the points, between T waveform and next period of P waveform, to 0. This 
process facilitates elimination of noise and adjusting the baseline of ECG signal to 0. 

3.0.3. Wavelet Transform Based Rerative Thresholding (WTIT) 

After preprocessing of ECG signal having completed through GBM, it is decomposed by discrete wavelet 
transformation (DWT), which involves passing this signal through high and low pass filters. The wavelet 
coefficients are obtained by the convolution process on the input signal with the impulse responses of high 
(g) and low (h) pass filters. Eqn.10 provides the output of high pass filter after performing downsampling on 
convolved signals to give detailed coefficients (CD), while Eqn.ll obtains the output of low pass filter that 
gives approximate coefficients ( CA ) after performing downsampling . 

OO 

CD[k] = x[n].g[2k — n] (10) 

n =—oo 

OO 

CA[k]= x[n].h[2k — n] (11) 

71 — — OO 

This WTIT algorithm is mainly dependent upon two system parameters: (i) Window Size ( Wi n ) and (ii) 
Threshold (T),,). How these parameters directly affect the performance of the energy aware sparse encoding 
scheme, discussed in the next few sections. A window is traversed along the detailed and approximate wavelet 
coefficients. The number of coefficients inside the window is referred to as the window size. If the largest 
coefficient inside a window is M, then all elements with value lesser than Th * M are approximated to zero. 
Here, Th is referred to as the normalized threshold. Its ranges vary from 0 to 1. In particular, we use WTIT 
to extract the high frequency components of ECG signal. This WTIT sets high approximated value for 
lower level coefficient, and lower approximation to higher level coefficient. This is continued to reach higher 
sparsity in ECG signal, keeping the higher frequency components intact. These coefficient matrices obtained 
from WTIT are transmitted after having undergone Huffman encoding. Huffman encoding is a lossless data 
compression technique which encodes sparse matrices in terms of only non-zero values. In the receiver end, 
the received signal has undergone Huffman decoding and inverse wavelet transformation for time domain 
representation of the signal to be analyzed by the person under concern. 
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4. System Design 

4-.1. Interaction between System. Parameters 

This energy efficient IoT-health monitoring system requires a mathematical trade-off between quality 
of the reconstructed signal and the energy spent in the Processing and Communication subsystems that is 
shown in Fig.l. As wavelet coefficients provide a sparse representation of a signal, wavelet based thresholding 
is applied on the ECG signal that is preprocessed by applying GBM. Furthermore, WTIT provides two 
independent parameters, threshold (Th) and window size (Wj n ). Th and Wi„ are used as approximation 
knobs for the entire system. The selection of sparse wavelet coefficients for representation of a signal is 
supervised by Th- Moreover, Wj„ can be adjusted according to the need of study of ECG transients for 
clinical purposes. So, Th and Wi n are the independent approximation parameters which control the IoT- 
based health monitoring system. 

Fig- ID and Fig. [7] have been drawn from experimental results. They represent the relationships between 
error bound and system energy consumption to the threshold parameter and the window size. When the error 
bound e r is the deviation between original and reconstructed signals, Fig. [6] represents how e r varies with Th 
for various Wi n of lengths 128, 256, 512, 1024, 2048 and 4096 respectively. If Th increases, the ECG signal 
encoded by WTIT becomes more sparse, and hence, e r increases. With an increase of Wj n , fewer coefficients 
are thresholded to zero, thus e r decreases. If mc energy is the energy consumed by the microcontroller during 
the operation and tx energy is the transmission energy then tot energy = mc energy + tx energy . and Fig. [ 7 ] 
represents the variation of tot energy with Th for various values of Wi n . To decrease the range of variation 
of tot energ y w.r.t Wi n and Th, logarithm scale is used. While mc en ergy and microcontroller execution time 
( mctime ) are nearly independent of Th, txenergy changes inversely related to Th- With an increase in the 
number of sparse coefficients, fewer bits are transmitted. Then, transmission time ( tx time ) decreases, which 
results in a decrease in tx ene rgy Thus, tot energ y decreases marginally with an increase in Th . However, as 
Wi n increases, the sparsity of encoded coefficients decreases and hence tx ener gy increases. 

Fig. [ 6 ] and Fig. [ 7 ] are obtained from experimental results that represent the nature of relationship between 
system parameters, thus driving us to develop an energy optimization heuristic. Therefore, provided a 
maximum tolerable error bound, we need to select that coordinate of window size and threshold parameter 
that result in the least total energy consumption. 

Table 2: Mathematical relations among various system parameters (c and a being positive real number variables) 


Parameters 

Curve fitting expression 
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Figure 6: Variation of \og 10 [tot energy ] with T) l and W ln obtained from experimental results 

4-1-1- Effect of window-size variation 

Upon variation of approximation parameters W and Th, the change of values for e r , mctime and txume 
are observed. These observed values are used to generate polynomial relations among the system parameters 
that are given in Table [2] These relations are generated by polynomial regression of the values of system 
parameters collected from actual experimentation executed in this work. The types of the polynomials clearly 
indicate that e r decreases as W m increases. This effect can also be observed from Fig. [6] When the length of 
window increases in WTIT, the less number of coefficients are thresholded to zero. This results the decrement 
of e r . The time complexity of WTIT is in the order of 0(n) and hence mctime increases linearly with the 
signal length under consideration. Increasing the Wi n results in improvement of quality of the reconstructed 
signal at the expense of increased mc energ y and tx energy . These are illustrated in both Fig. [6] and Fig. [7] and 
implied from Table [2] 

4-1.2. Effect of threshold variation 

To study the effect of variation of Th, e r is varied with Th for a fixed IIy rt ■ With increase in the value of 
Th, the number of sparse wavelet coefficients increases, thus introducing subtle differences between original 
and reconstructed signals. From Fig. [6j we see that e r varies linearly with Th- As Th increases, fewer 
coefficients are transmitted as rest of the coefficients below the threshold level are set to zero. These zero 
values, discarded via Huffman encoding, are not transmitted. The non-zero coefficients are transmitted in 
the form of packets. Higher the number of non-zero coefficients the larger is the packet-size. Thus, tx t imei 
and hence tx energy , are directly proportional to the packet-size and inversely proportional to Th- However, 
Th does not affect mctime and mc ener gy It is observed that Th has negligible influence when tot energv is 
calculated. In Fig. [7J tot energy has a small negative slope, when varied with Th for a given W; m because 
tx energ y decreases with the increase of Th- 

Thus, the IoT-based ECG monitoring system is designed in such a way that for a pre-specified maximum 
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Figure 7: Variation of e r with Th and Wi n obtained from experimental results. 


tolerable error bound ( K ), the tot energ y is to be minimized by varying the independent approximation 
parameters Th and Wi n . This requirement of IoT-based ECG monitoring system demands to solve an 
optimization problem. 

4-2. Optimization Heuristic 


Table 3: Notation of system parameters 


Notation 

Sub-parameter variable 

W in 

T h 

txtime 

fqr 

energy 

'TICtime 

TlC ener gy 

t , Ot erier gy 

Gr 

window-size 
threshold parameter 
module transmission time (ms) 
module transmission energy 
microcontroller system time 
module microcontroller energy 
total energy ( tx energv + mc energy ) 
error between original and reconstructed signal 


Table 4: Sum of squared errors (SSE) using thin plate spline interpolant 


Modelled function 

Interpolant type 

Goodness of fit 

e r (W in ,T h ) 
t , Ot e nergy{y^in-) Th) 

Thin plate spline 
Thin plate spline 

SSE = 2.78e-27 
SSE = 1.625e-08 


In order to discuss the optimization heuristic, the notations for sub-parameter variables are given in Table 
3. Variables e r , tot energy , mc energy , mc t i me , tx energy , tx t ime are dependent on the independent variables W in 
and Th- Both the variables e r and tot energy can be represented as a function of independent variables Wi n 
and Th . From simulation data, the observed values are modelled to form 3-dimensional curves by using 
curve-fitting toolbox in MATLAB 2017a. Thin plate spline is chosen as the interpolant because it gives the 
best fitting model with least sum of squared errors (SSE) as compared to other best fitting models like linear, 
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polynomial and power fitting. The ideal fitting will have a SSE tending to zero. The goodness of fit using 
thin plate spline is given in Table [4j 

To develop an energy efficient system, we minimize the tot energy for a given maximum tolerable error 
bound {K). So, the optimization problem can be stated as : 

Minimize tot energ y(W in ,T h )-, 
subject to e r ( W ln ,T h ) < I\ 

The error bound K is set as preset requirement for studying transients of ECG signal. The optimization 
heuristic namely Algorithm 1 chooses the optimum ( Wi nopt ,Th opt ) coordinate from the set of coordinates 
(W lnt , Tfr. ) where i £ Z. Algorithm-1 is a optimum parameter search algorithm that does not take into account 
the parametric relationship between the independent variables ( Wi n and XK) and the system parameters (e r 
and tot energy ). In Fig-10, the variation of Total Energy with Window-size is much higher than the variation 
of Total Energy with Threshold Parameter. Exploiting this fact Algorithm-2 is developed which reduces the 
algorithm complexity from 0(n 2 ) to 0(n). 

J^.2.1. Selection of coordinates that satisfy inequality constraint 

The set of coordinates (W lrli . Tf, . ) that satisfy the inequality constraint of the optimization problem is 
denoted by set G. All the operating points which belong to the set G has value of Normalized Error less than 
the error bound K. Hence, (Hy rti ,T) ti ) £ G implies ( Wi Ui ,Th i ) < K. Selection of operating points in set G 
will be shown in Fig. 9 in Section 6. 

f.2.2. Selection of optimum coordinate from set G 

The search space for the optimum coordinate is narrowed down by selecting operational points from set 
G. The coordinate fWi ni ,Thf) that minimizes the value of tot energy is the optimum coordinate. Algorithm 
1 applies brute search method on set G to obtain the optimum operating point (W ,;„ opt ,). Selection of 
the optimum operating point from the set G will be shown in Fig. 10 in Section 6. 

4-2.3. Deriving the Quality-Energy (QE) relation 

The previous two steps are repeated for various values of error bound K, and consequently, the optimum 
coordinate is obtained for each value of K. This optimum coordinate minimizes total energy and is chosen for 
energy-efficient operation of IoT based health monitor system. Whenever this IoT based health monitoring 
system is in operation for some given values of (K). the system operates at (Wi Uopt , Th opt ), which is the optimal 
operating point for the minimization of the total energy. In this case, Fig. [TT] will show the Quality-Energy 
(Q-E) relationship of the system that denotes the minimum achievable energy for a given error bound. 

Algorithm 1 searches for all the coordinate points in set G making it of complexity Oin 2 ). However, we 
have already obtained a relation of Total Energy with both Threshold Parameter and Window-Size in Fig. 
10 . Exploiting this fact that variation of tot energy with respect to Wi n is greater than, we first perform a 
linear search on window size and obtain the Wi nopt . Then keeping Wi Uopt fixed, we perform a second linear 
search on the Threshold Parameter to obtain Th opt . This reduces/prunes the search space significantly as 
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Algorithm 1: Finding optimized values of Wi n and Th for energy minimization (Wi nopt , Th opt ) using a 
generic search method 


Input: Error Bound K 

Output: Optimum window-size and threshold (Wi nopt ,Th opt ) 

1 G = 0 ; 

2 for (i = 0;i< n; i=i+l ) do 


// G is a null set 


while e r (W ini ,T hi ) < K) do 


G =G U 


// include the new coordinate into the existing set 


end 


4 

5 

6 end 

7 g = |G|; // cardinality of the set G 

8 

9 for (i = 0; i< g; i=i+l) do 

10 (W ini ,T hi ) = (oo,oo) 

11 if tot ener gy {XVi ni ? Thi ) — totenergy (w inop , t > T hopt ) , 

12 (Win 0 t ,T ho t ) = (W ini ,T hi ) ; // select the optimum coordinates from the set G 

13 end 

14 return (W iriopt , T hopt ) 


Algorithm 2 reduces to a 0(n) search problem as when compared to Algorithm 1 which is a 0(n 2 ) search 


problem. Hence, by exploiting the relationship of independent parameters with Total Energy, we avoid a 


generic search algorithm. This allows us to extract the final design parameters of Wi nopt and Th opt quickly 
and are the same as the ones shown (emboldened) in Table [5] 


5. Experimental Methodology 

In order to describe the hardware experimental set-up, we refer to Fig|4j We use AD8232 ECG heart 
rate monitor with an Arduino UNO R3 for ECG signal acquisition from a subject (patient) body. The set of 
pins of AD8232 are connected to the corresponding pins of Arduino; those pins are namely ground pin, 3.3V 
power supply pin, output signal pin, positive leads-off detect pin, negative leads-off detect pin of AD8232 
which are connected with ground pin, 3.3V pin, AO pin, pin number 10 and pin number 11 of the Arduino 
respectively. AD8232 has three coloured cables: black, blue and red, which are placed on the right arm, left 
arm and right leg of the patient respectively. Arduino C processing IDE v2.2.1 is used to write sketch to read 
and process ECG signals directly. Arduino IDE shows a list of available communication ports, any of which 
can be used for ECG signal recording. From the Arduino, ECG signal is transmitted to a server located in 
private/public cloud using ESP826 module which is used to access the Internet via WiFi connection. The 
energy consumption of the system is measured using a fluke high precision multimeter. 

We apply static thresholding, i.e., it is assumed that during the transmission period, the heart condition 
of the patient remains unchanged or the threshold level is constant. However, the physician can change 
the normalized error value on ECG monitor as and when it is required to study the ’Quality’ of the signal 
(refer to Fig. 11). We perform the experiment by taking a generalised value of Th, which will be further 
changed. In our simulation, MATLAB R2017a is used to evaluate the mathematical relations among various 
system parameters that are given in Table 2. We use the Curve Fitting Toolbox in Matlab R2017a, which 
would interpolate the sample observation points to generate a well-defined function in order to study the 
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Algorithm 2: Finding optimized values of W,; n and Th for energy minimization ( Wi Uopt , Th opt ) exploiting 
parametric relationship 


Input: Error Bound K 

Output: Optimum window-size and threshold (W in ovt ,Th ovt ) 


1 G = 0 ; 

2 for (i = 0;i< n; i=i+l ) do 

3 while e r (Wi ni ,Tf li ) < K) do 

4 | G = Gu'(W i n i ,T hi ); 

5 end 

6 end 

7 9 = |G| ; 

8 for (ii = 0; i< (/; i=i+l) do 

9 
10 


// G is a null set 

// include the new coordinate into the existing set 

// cardinality of the set G 


C Win opt ) = (oo); 

if tot ener gyfWi ni ) <:tot ener gy (Win opt ') , 

Wi nopt =Wi ni ; // select the optimum window-size coordinate from the set G 

12 end 

13 for (i = 0 ; i< q; i=i+l) do 

14 ( T K pt ) = (oo); 

15 if tot ener gy (y^in opt 5 Thi ) <tot ener gy (WiriopuThopt) 1 

16 Th opt =Th i ; II select the optimum Threshold Parameter coordinate 

17 end 

is return (W inopt ,T hopt ) 

dependencies among various system parameters. The interpolating function is chosen in such a way that it 
minimizes SSE of the observed points and the interpolating curve. 


6. Experimental Results 

In this work, we have developed a real time encoding scheme which performs iterative thresholding and 
approximation of wavelet coefficients for sparse encoding of ECG signals. We have optimized various process 
parameters to model a low power hardware prototype for the implementation of our algorithm on a real time 
microcontroller based IoT platform that exhibits an end-to-end WBSN system. 

In order to describe the performances of these algorithms, various system parameters and finally the 
implementation of hardware prototype, we start the discussion the results from Fig. [8j The nature of 
variations of rnctime with W-i n is evident in Figjsji. With the increment of W ln , mctime increases almost 
linearly. This graph depicts that when Wi n increases from sample size window of 16 to sample size window of 
4096, mctime increases from 0 ms to 73 ms. This is because of the fact that the asymptotic time complexity 
of WTIT algorithm is in the order of O(n). As lF,; n increases, the more number of discrete samples of ECG 
signals under consideration are processed together by WTIT, and hence mctime increases. However, mctime 
remains constant with respect to T^. Irrespective of the variations of Tf ll the number of discrete samples 
processed by WTIT remain unchanged. Hence, mctime does not vary with the variations of Th- 

Fig. H, and Fig. describe the variations of e r with sub-optimal values of Wi n and Th respectively. For 
a fixed value of Th, as Wi n increases, e r decreases due to the fact that the WTIT algorithm approximates 
fewer wavelet coefficients to zero. When W ln increases from 512 to 4096 sample sizes, e r decreases from 0.043 
to 0.02. An increase in Th results in higher number of sparse wavelet coefficients. Normalized e r increases 
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(a) mc time vs W in 



(b) 


(b) e r vs W in 



W,„ 


W 


( c ) (c) e r vs T h 



T h 


Figure 8: Dependencies of various subparameters on independent metrics obtained from experimental results : (a) Variation of 
microcontroller execution time and window size (6) Variation of normalised error with window size (c) Variation of normalised 
error with threshold 

from 0.1 to 0.15, when Th increases from 0.01 to 0.015. Thus, e r and T\ are directly proportional to each 
other. Similar nature of variations of e r with Th and Wi n was illustrated in FigjGj 

Fig. i demonstrates the relationship of e r with Wt n and Th using a 3-D surface generated by thin plate 
spline interpolation having used on the experimental data. On the variation of both Wi n from 16 to 4096 
sample sizes, and Th from 0.005 to 0.02, normalized e r varies from 0.05 to 0.2. To explain Fig. [lUJ it is 
imperative to discuss that txume increases and decreases linearly with respect to Wi n and Th , respectively. 
With the increment of sample sizes in Wi n , WTIT has processed these number of samples of ECG signal 
and subsequently transmitted to server located in private/public cloud.This results in the increase of txume- 
As Th increases, the number of sparse wavelet coefficients increases. Thus, after having undergone Huffman 
encoding, both the amount of transmitted data and txume decrease. Note that txume and rriCtime are directly 


related to tx 


energy 


and mc ener gy respectively. As defined earlier, tot energy is the sum of mc energy and tx energy . 


It is observed that while mc er 
results in the negative slope of tot, 


does not vary with Th , tx e 


decreases with an increase of Th . This 


energy versus Th that was illustrated in Fig. [T] Since both mc energy 


and 


tx 


energy 


increases with W ln , tot, 


energy 


increases with Wi n . In Fig. [Tol generated with experimental data, 


similar variation of tot„ 


with different values of T h and W in is observed in the 3-D surface generated by 


thin plate spline interpolant. With the variation of both Wi n from 16 to 4096 window sample sizes, and Th 
from 0.005 to 0.02, tot P 


^energy 


varies from 0.008 mJ to 4.71 rnJ. It is evident in Fig. 10 
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Figure 9: Variation of normalized error with window size and threshold obtained from experimental results. A red plane 
indicating constant error of K is drawn which cuts the surface. All coordinates lying below this red plane belong to the set G. 
Algorithm 1 and Algorithm 2 search for (W ln<lpl , T hapl ) for all values of (W irl .. T hj ). From Fig. |9j given an 

error bound I\ marked by the red plane, the points (W lrii , Tt H ) are chosen in such a way that they lie below 


this plane. For a given (Wi nj , T/jJ pair, tot energy is determined that is shown in Fig. 10 The pair ( Wj ;n% , T)j.) 
for which tot energy becomes minimum, is the optimal operational point (W,:„ opt , Th t ) for the design of IoT 
based monitoring system. 

Table [5] enlists all {Win, , Th . L ) for various values of normalized error bound K. For each (W rni , T) ti| ), 
tot energy is calculated from the sum of mc en ergy and txenergy The optimal operation point ( Wi n . opt , Th opt ) 
minimizes tot ene rgy (marked boldfaced in Table [BJ for each K. For a normalized error bound A'=0.03, 
( w in opt , T h opt ) = (256,0.053) and the tot energ y consumed is 0.014 mJ. Similarly, for K = 0.05, {W iriopt ,T hgpt ) = 
(128, 0.0072) and the tot energy consumed is 0.0101 mJ. For various values of K (marked boldfaced in Table 
[5]), the optimal operational point {Wi Hopt , Th opt ) is derived for all values of (W. in .. ) using Algorithm 1 and 

Algorithm 2. Accordingly, the tot energy consumed is calculated and given in Table [5j Fig. [TT| shows the 
final Quality-Energy (QE) relationship of quality and energy parameters and it demonstrates a significant 
decrease in optimized tot en ergy consumption for a negligible degradation of the quality of signal. Through 
the application of Algorithm 1 and Algorithm 2, energy consumption is reduced to 96% with a deterioration 
of 2% in signal quality. Also when the deterioration of the signal quality is 1%, the energy consumption is 
reduced to only 22% as shown in Fig. 11. This abrupt change in the energy consumption due to change in 
error level leads to the selection of the optimum Normalized Error bound by elbow method shown in Fig. 
11.Thus, energy consumption decreases with a increase in error bound K. With an increment in tolerable 
error, sparsity of the encoded vector increases resulting in a significant amount of the decrease of tx en ergy 

We introduce the quality/error-levels (degree of tolerable error) just for demonstrating the benefits of 
approximation of the ECG signal. These quality/error levels can be adjusted by the physicians as per 
the quality of ECG signal originated from patients and then these values can be adapted accordingly. To 
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Figure 10: Variation of total energy with threshold and window size obtained from experimental results. The coordinates 
belonging to the set G are shown by the arrowed region. The optimum coordinate is the one that simultaneously belongs to 
this region and has the least value of total energy. 

demonstrate that Energy-Efficient compression-transmission scheme can be generalised to any standard ECG 
signal, we derived the Q-E plot for three different ECG signals namely ECG-1, ECG-2 and ECG-3 respectively 
obtained from three different patients each having a different health condition. This leads to the conclusion 
that the irrespective of the health conditions of the patients, the compression-transmission scheme for those 
ECG signals shows approximately 96% energy savings with negligible deterioration in signal quality (2%). 

7. Conclusion 

In this work, an energy efficient hardware prototype developed on a IoT-based microcontroller platform 
has been deployed to operate as an end-to-end WBSN system in real time. This WBSN has the capability 
of the compression of real time ECG signal via sparse encoding of these signals in time-frequency domain 
using wavelet coefficients. The signal compression scheme has been developed a two step approach: GBM has 
encoded the signal in time domain and WTIT has worked in time-frequency domain to perform approximation 
of wavelet coefficients. However, excessive approximations of these coefficients result into an appreciable 
degradation of signal quality and the loss of clinically important information. On the other hand, sub- 
optimal amount of approximations has generated large volumes of transmission data that consumes higher 
energy and higher bandwidth. In this work, we have proposed a methodology to perform optimum trade-off 
between the quality of signal, the consumption of energy and the usage of bandwidth (shown in Fig. 11). 
Algorithm 1 and Algorithm 2 have optimized the various design parameters by introducing optimal amount 
of approximations that have led to develop low power IoT-based ECG health monitoring hardware prototype. 
With these proposed algorithms, energy consumption of the IoT system has been reduced by 96% with only 
2 % deterioration in signal quality. 
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Figure 11: Q-E relationship diagram denoting the minimum achievable normalized energy for a pre-specified error bound K. 
Normalized energy is Total Energy consumed in our compression-transmission scheme divided by the Total Energy without any 
prior compression. Using the Elbow method, the optimum Normalized Energy is found to be 0.04 indicating decrease in energy 
consumption of about 96%. Q-E plot has been obtained for three different ECG subjects 


Table 5: Variation of optimum system variables for different values of K 


K 

(W ini ,T hi ) 

mctimei ms) 

TnC-energy (mj) 

txume (ms) 

energy (nij) 

t0t er i er gy (mj) 

0.015 

( 2048 , 0 . 0056 ) 

191.98 

0.37 

55.07 

0.2297 

0.6069 


(4096,0.0071) 

823.60 

1.61 

742.23 

3.096 

4.71 

0.02 

( 1024 , 0 . 0054 ) 

41.71 

0.08 

3.21 

0.0134 

0.0953 


(2048,0.0076) 

191.99 

0.3772 

54.99 

0.2294 

0.6066 


(4096,0.0096) 

823.66 

1.618 

741.26 

3.092 

4.71 

0.03 

( 256 , 0 . 0053 ) 

1.42 

0.0028 

2.68 

0.011 

0.014 


(512,0.0134) 

7.57 

0.0156 

2.61 

0.0109 

0.0265 


(1024,0.0176) 

41.71 

0.082 

3.20 

0.0134 

0.0954 


(2048,0.0116) 

192.01 

0.3772 

54.83 

0.228 

0.605 


(4096,0.0149) 

823.7 

1.61 

739.21 

3.08 

4.70 

0.04 

( 128 , 0 . 0057 ) 

0.745 

0.00146 

2.069 

0.00863 

0.0100 


(256,0.0071) 

1.424 

0.0028 

2.679 

0.0117 

0.0145 


(512,0.0087) 

7.95 

0.0156 

2.605 

0.0108 

0.026 


(1024,0.0114) 

41.71 

0.081 

3.195 

0.0133 

0.095 


(2048,0.001575) 

192.03 

0.377 

54.67 

0.228 

0.605 

0.045 

( 128 , 0 . 0064 ) 

0.746 

0.00147 

2.065 

0.00861 

0.01008 


(256,0.0089) 

1.42 

0.0028 

2.67 

0.011 

0.0139 


(512,0.0098) 

7.95 

0.00156 

2.60 

0.010 

0.0264 


(1024,0.0129) 

41.71 

0.081 

3.19 

0.0133 

0.0952 


(2048,0.018) 

192.04 

0.377 

54.58 

0.227 

0.6055 

0.05 

( 128 , 0 . 0072 ) 

0.745 

0.0015 

2.069 

0.0086 

0.0101 


(256,0.0089) 

1.42 

0.0029 

2.67 

0.011 

0.014 


(512,0.010) 

7.95 

0.0156 

2.60 

0.010 

0.0264 


(1024,0.0144) 

41.72 

0.082 

3.18 

0.013 

0.095 

0.06 

( 64 , 0 . 0059 ) 

1.05 

0.0021 

1.50 

0.0062 

0.008 


(128,0.0087) 

0.7449 

0.0015 

2.05 

0.0085 

0.010 


(256,0.010) 

1.42 

0.0028 

2.66 

0.011 

0.013 


(512,0.0134) 

4.95 

0.0097 

2.58 

0.0108 

0.0205 


(1024,0.0176) 

41.722 

0.082 

3.17 

0.0132 

0.095 
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